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ABSTRACT 

We present spectroscopic confirmation of tlie CI 1604 supercluster at z ~ 0.9. Originally detected 
as two individual clusters, CI 1604+4304 at z = 0.90 and CI 1604+4321 at z = 0.92, which are 
closely separated in both redshift and sky position, subsequent imaging revealed a complex of red 
galaxies bridging the two clusters, suggesting that the region contained a large scale structure. We 
have carried out extensive multi-object spectroscopy, which, combined with previous measurements, 
provides ^ 600 redshifts in this area, including 230 confirmed supercluster members. We detect two 
additional clusters that are part of this structure, CI 1604+4314 at z = 0.87 and CI 1604+4316 at 
z = 0.94. All four have properties typical of local clusters, with line-of-sight velocity dispersions 
between 489 and 962 km s~^. The structure is significantly extended in redshift space, which, if 
interpreted as a true elongation in real space, implies a depth of 93 h^Q Mpc. We examine the spatial 
and redshift distribution of the supercluster members. 

Subject headings: catalogues - surveys - galaxies: clusters: general - large-scale structure of the 
Universe 
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1. INTRODUCTION 

The CI 1604 supercluster was initially detected as 
two separate clusters, CI 1604+4304 at z = 0.90 and 
CI 1604+4321 at z = 0.92, in the plate-based survey 
of iCunn. Hoessel. & Okj ((l986). Deeper imaging and 
multi-object spectroscopy taken with the Low Resolu- 
tion Imaging Spectrogr aph (LRIS; Okc ct al . 1995) at the 
Keck telescopes by Ok e. Postman, fc LubinI |1998, here- 
after 098) yielded redshif ts and preliminary velocity dis- 
persi ons for each cluster ijPostman. Lubin. fc 
1211011, hereafter P98,P01). Motivated by the close sepa- 
ration of CI 1604+4304 and CI 1604+4321 in both radial 
velocity (43 km s""'^) and position on the sky (17'), 
iLubin et al.l l|2000ri performed deep multi-band imaging 
with the Palomar 5-m telescope covering the area be- 
tween the two clusters. The imaging revealed an over- 
density of red galaxies whose colors were consistent with 
early-type galaxies at z ~ 0.9, suggesting that the clus- 
ters were part of a high-redshift supercluster. 

To verify this conclusion, map the supercluster struc- 
ture, and obtain a large sample of cluster galaxies with 
measured spectroscopic properties, we are conducting an 
extensive spectroscopic survey spanning the region be- 
tween CI 1604+4304 and CI 1604+4321. In this Let- 
ter, we present spectroscopic confirmation of this large 
scale structure and a discussion of its properties based 
on the currently known supercluster members. We as- 
sume a ACDM cosmology with flm = 0.3, A = 0.7, and 
Ho = 70 hm km s'^ Mpc"!. 

2. THE SPECTROSCOPIC SURVEY 

Here we briefly describe the spectroscopic survey of 
the CI 1604 supercluster. A full discussion of the data 
acquisition and analysis procedures will be presented in 
the catalog paper (Gal, Lubin & Oke 2004, in prepa- 
ration). The survey was conducted using LRIS and 
the Deep Imaging Multi-object Spectrograph (DEIMOS; 
iFaber et al.ll2003j ) on the Keck 10-m telesc opes. Target 
selection was based on multi-band imaging ijLubin et al.l 
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|2000() taken at the Palomar 5-m telescope. A contiguous 
area of 10.4' x 18.2' (4.8 Mpc x8.5 Mpc) between the 
original two clusters was imaged in _B, R, and Gunn i 
filters. The galaxy catalog reaches a depth of i? ~ 26, 
i? ~ 25 and z ~ 24 for a 5cr detection. The resulting 
color-magnitude diagrams reveal a well-delineated red se- 
quence of galaxies whose colors (1.2 < R — i < 1.7) are 
consistent with spectroscopically-confirmed, early-type 
galaxies in the tw o original clusters (see Figures 2 & 4 of 
ILubin et alJl200(Tl) . In Figure we indicate on the com- 
posite i band image all of these galaxies by overlaying 
contours of red galaxy density. These contours clearly 
delineate the large system of galaxies encompassing the 
original clusters. We also observe two strong concentra- 
tions of red galaxies directly between CI 1604+4304 and 
CI 1604+4321, suggesting the presence of two additional 
clusters. Objects within the red sequence were preferen- 
tially (but not exclusively) targeted for slit assignment in 
our current spectroscopic survey. Targets were observed 
with a combination of LRIS and DEIMOS slitmasks. 

LRIS slitmasks covering parts of this region were ob- 
served on UT May 4-7, 2000. A total of six slitmasks 
were used, with ~ 30 objects per mask. Objects with 
magnitudes down to i — 23.0 were included. Data were 
taken using the 400 1/mm grating blazed at 8500 A, with 
a central wavelength of 7000A. This setup provides spec- 
tral coverage from 5500A to 9000A and a dispersion of 
1.86A pix ^; with the 1" slits utilized, the resolution is 
7.8A (550 km s'^ at z = 0.9). Seeing was typically 0.7". 
Four of the masks were exposed for 4 x 1200s each, one 
mask for 3 x 1200s, and one for 2 x 1200s. HgNeArKr 
arc lamps, as well as flat fields, were observed immedi- 
ately before and/or after each set of target exposures to 
provide wavelength calibration with the instrument and 
telescope in the same configuration, minimizing errors 
due to flexure. The data were processed using standard 
IRAF tasks and scripts. 

To obtain a larger sample of objects, extending to 
fainter galaxies, and covering a wider spatial area, further 
spectroscopy was obtained with DEIMOS on UT 25-26 
May 2003. Galaxies as faint as z = 24 were selected, as 
well as some objects previously observed with LRIS but 
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Fig. 1.— Composite 10.4' x 18.2' (4.8 Mpc x8.5 Mpc) i band 
image of the supercluster (N is up, E is to the right) Overlaid are 
contours of red (1.2 < R — i < 1.7) galaxy density. The original 
clusters CI 1604+4304 (A) and CI 1604+4321 (D) are at the bottom 
and top of this image. The new clusters CI 1604+4314 (B) and CI 
1604+4316 (C) are detected as significant overdensities. Galaxy 
density contours range from to 9 galaxies arcmin"^, in steps of 
0.75 galaxies arcmin"^. 



for which no redshift was obtained. We targeted ^ 85 
objects on each mask. The instrument was configured 
with the f 200 1/mm grating, blazed at 7500A. The cen- 
tral wavelength was set to 7700A such that OIIA3727 and 
Ca H and K AA3934, 3968 at the redshift of the super- 
cluster would be easily visible. Slit widths were 1", and 
where possible, slits were aligned with the major axis of 
the galaxy. This configuration results in a pixel scale 
of 0.33A pix ^, a resolution of '^1.7A (120 km s ^ in 
the cluster frame), and spectral coverage from 6385A to 
9OI5A. We observed four masks with a total exposure 
time of 12600s (7 x 1800s) per mask. The large num- 
ber of exposures permits excellent cosmic ray subtraction 
and a significant reduction in background noise. Seeing 
ranged from 0.9" to 1.5" on the first night and was nearly 
constant at ~ 0.6" on the second night. For each mask 
setup, we obtained a set of three internal flatfields, as 
well as an arc line calibration exposure with the NeAr- 
KrXe la mps. Data reduct ion was performed using the 
DEEP2 (Davis et alJ 12003(1 version of the spec2d data 
reduction pipeline. Based upon the Sloan Digital Sky 
Survey (SDSS) spectral reduction package, this IDL soft- 
ware performs all of the standard processing steps, pro- 
ducing sky-subtracted, wavelength-calibrated one- and 
two-dimensional spectra^. 

Redshifts for the LRIS and DEIMOS observations were 
all measured using an IRAF task called redsplot (T. 
Small, private communication). This tool allows the user 

^ For det ails on this pipeline, see the DEEP2 website at 
|http : //astron . berkeley . edu/~cooper/deep/spec2d/ 
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Fig. 2. — The redshift distribution for all galaxies within 0.8 < 
2 < 1.0. The four clear peaks correspond to the four member 
clusters in the supercluster. Each component cluster is shaded and 
labeled with its identifying letter (see Table [Tl . 



to interactively identify features in the one-dimensional 
spectra and displays the locations of other possible fea- 
tures from a user- provided spectral line list. In this way, 
one can centroid various visible features and test differ- 
ent redshifts to find which produces a consistent set of 
identified lines. Median velocity errors for the LRIS and 
DEIMOS data are 150 and 90 km s^^, respectively. 

The six LRIS masks contained a total of 179 slitlets; 
of these, 156 were distinct targets. An additional four 
objects were found serendipitously for a total of 160 
unique objects. Identifications were made for 147 of 
these, with 116 (80%) being extragalactic, and 65 of these 
(56%) having redshifts within the supercluster range of 
0.85 < z < 0.96. The four DEIMOS slitmasks had 326 
unique slitlets. An additional 57 objects were detected 
serendipitously, yielding a total of 383 objects. Of these, 
we were able to identify 309 objects, only 22 of which are 
stars. Over 40% (122) of the identified objects were part 
of the supercluster complex. The combined LRIS and 
DEIMOS data give 378 unique extragalactic redshifts. 

Our analysis also includes the original sample of 
LRIS spectroscopy measured for CI 1604-1-4304 and CI 
1604+4321 (098). Objects observed as part of that sur- 
vey were distributed over an approximately 2' x 7' region 
centered on each cluster. Redshifts were measured for 
103 and 135 galaxies in the fields of CI 1604-1-4304 and 
CI 1604+4321, respectively (P98,P01). 

The combined spectroscopic survey yields a total of 
597 unique extragalactic redshifts. In Figure |21 we show 
the redshift histogram for the 289 galaxies in the range 
0.8 < z < 1.0. There are four clear peaks in this dis- 
tribution which correspond to the four member clusters 
of the supercluster. Redshift peaks at z = 0.90 and 
z = 0.92 belong to the original clusters, CI 1604-1-4304 
and CI 1604-f4321, respectively. The new clusters iden- 
tified through the multi-band imaging (see Figure^ cor- 
respond to the redshift peaks at z = 0.87 (CI 1604-f 4314) 
and z = 0.94 (CI 1604-1-4316). Table|T]lists the name, the 
reference letter used throughout the Letter, the central 
coordinates, and the median redshift for each cluster. 
The central coordinates and median redshift are mea- 
sured from those galaxies which are counted as members 
in the detailed analysis of cluster membership presented 
in §3. 

Figure |31 shows the location of all 597 spectroscopic 
galaxies, with the 229 confirmed supercluster members 
coded by cluster membership. Although the galaxies be- 
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Fig. 3. — The spatial distribution of 597 galaxies with redshifts in 
the CI 1604 supercluster region, using all of the LRIS and DEIMOS 
data. Supercluster members are plotted with large symbols, coded 
by cluster membership, while non-members are plotted as points. 
Dashed circles represent a 1 h^^ Mpc (projected) radius centered 
on each cluster. 



longing to each member cluster appear reasonably con- 
centrated and separated from one another, there is sig- 
nificant overlap throughout the supercluster region. In 
addition, the spatial distribution of galaxies suggests that 
we have not yet reached the edge of this structure. We 
note, however, that the region has not been uniformly 
sampled. 

3. CLUSTER PROPERTIES 
3.1. Velocity Dispersions 

Velocity dispersions are notoriously difficult to calcu- 
late for galaxy clusters, although this has not prevented 
the publication of results using as few as ten redshifts 
per cluster. Selecting cluster members, a radius in which 
to include objects, the mathematical technique, and even 
the mean cluster redshift each pose their own subtle is- 
sues. These problems are compounded in densely pop- 
ulated areas where the individual clusters overlap in a 
larger structure. Nevertheless, we present initial esti- 
mates of the cluster velocity dispersions, while postpon- 
ing detailed discussion for a future work, with the goal 
of demonstrating that the supercluster constituents are 
indeed significant structures, as suggested by POl. 

To determine the membership of each superclus- 

ter component, we utili ze the ROSTAT package 

(|Beers. Flvnn. fc Gebhardl IT990|) . which measures a 
wide variety of estimators for the location (median) and 
scale (line-of-sight dispersion) of a set of measurements. 
It avoids assumptions of gaussianity in the underlying 
velocity distributions (which is unlikely to be true for in- 
teracting clusters) and includes bootstrap and jackknife 
error estimation for the cluster redshifts and dispersions. 
Many of these estimators are significantly more resistant 



to non-gaussianity in the sample than the traditional 
Gaussian; large differences between the different scale 
estimates can be a sign of deviations from gaussianity 
(either in shape or due to outliers) in the sample. 

We first visually divide our redshift histograms into 
separate (but overlapping) redshift ranges for each clus- 
ter, selecting a conservative (broad) range for each which 
clearly includes the entire structure. The median redshift 
of each cluster is measured, and initial, cosmologically 
corrected peculiar velocities for each galaxy in the clus- 
ter are calculated. These velocities are input into ROS- 
TAT, which computes the biweight median velocity, and 
the line-of-sight dispersions abiwt and a gauss , among oth- 
ers. We then apply a shift to the velocities of each galaxy 
such that the biweight median velocity is zero, while pre- 
serving the dispersion. Objects with \v\ > 'iauwt are 
then clipped from the sample, and the median and dis- 
persions recalculated. This procedure is iterated until 
either the various dispersion measures converge, or there 
are no more objects to be clipped. For all clusters, the 
clipping results in removal of fewer than 10% of the in- 
put galaxies. The shaded regions in Figure |2 show the 
final redshift range for each cluster. Inspection of the 
spatial distribution of the clipped galaxies reveals that 
they are typically at large clustocentric distances. We 
note that using only those galaxies within the central 1.5 
Mpc yields comparable dispersions. 

For each cluster. Tabled lists the final redshift range, 
median redshift, and number of galaxies, as well as 
two estimates of the linc-of-sight velocity dispersion (in 
km s~^): Ubiwt, the biweight estimator and its 1-sigma 
bootstrap error l|Beers. Flvnn. fc Gebhar dt 1990), and 
gauss, the Gaussian estimator and its 1-sigma Student's 
t error. We note that the bootstrap error estimates are 
typically 25% to 40% higher than those estimated using 
the Student's t distributi on, which assumes an under- 
lying Gaussian population (jPanese. de Zotti. fc di Tullid 
|1980|1 . The aios estimates for clusters A and D are ~ 30% 
lower than those measured by POl, who used a factor of 
> 2 fewer galaxies. Comparison of our redshift histogram 
with their Figure 3 shows that POl included members 
of the distinct structures B and C when measuring the 
dispersions for A and D, respectively. Such discrepan- 
cies highlight the fact that velocity dispersions of clusters 
(and their error estimates) based on small (~ 20) redshift 
samples are subject to strong systematics and should be 
regarded with skepticism, especially in the presence of 
larger scale structures. 

3.2. Separations 

One goal of this study is to determine whether or not 
these clusters are components of a larger system (a super- 
cluster) or a chance projection on the sky. We estimate 
the cluster centers (listed in Table by starting with 
the position of the peak in red galaxy density for each 
cluster (see Figure QJ. We then calculate the median po- 
sition of all spectroscopically confirmed members within 
a radius of 1 h^^ Mpc. Under the assumption that these 
cluster are no longer expanding with the Hubble flow, 
we measure separations, based on these cluster centers, 
using the angular diameter distance. The largest trans- 
verse separation is 7.62 h^^^ Mpc (cluster pair AD), with 
the closest pair separated by only 1.85 h^^ Mpc (BC) 
in projection. The elongation in redshift space implies 
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Table 1. Cluster Properties 



Cluster 


ID 


RA 


Dec 


Rcdshift 


z range 


N 1 

gals 


'^los ,biwt 


^los, gauss 


CI 1604+4304 
CI 1604+4314 
CI 1604+4316 
CI 1604+4321 


A 
B 
C 
D 


16 04 23.7 
16 04 25.8 
16 04 05.7 
16 04 35.3 


+43 04 51.9 
+43 14 23.0 
+43 15 52.8 
+43 21 01.6 


0.9001 
0.8652 
0.9350 
0.9212 


0.8852 - 0.9146 
0.8539 - 0.8791 
0.9302 - 0.9420 
0.9111 - 0.9302 


67 
62 
24 
77 


962 ± 141 
719 + 71 
489 ± 79 
640 ± 71 


989+^^ 
726+11 

6491^9 



much greater radial distances, ranging from 18.1 (CD) 
to 93.0 (BC) Mpc. 

If we take the largest radial separation as the maximum 
extent of the supercluster, its size is not unusual for local 
superclusters (Bahcall 1988) and is consistent with the 
cluster-cluster correlation function, w hich exhibits power 
on scales as large as 100 /t~^Mp c (lBa.hcall fc Soneira 



on sc ales as large as lUU n Mpc aijahcall fc honeirg 
198,11 iPostman. Huchra. fc CTelieil 119921: iBahcall et al 
20031) . The large aspect ratio suggests 



that we are 

observing a filamentary structure nearly face-on, with 
the redshift differences too large to be due solely to 
peculiar velocities. Such elongated structures, with 
typical scales of ~ 50 h^Q Mpc, and as large as 100 h^Q 
Mpc, arc seen in large scale structure simulations 
j^atbv aprakash . Sahni. & Shandarin 1998; Evrard et al. 
20021 iKolokotronis. Ba silakos. & Plionis 2002), large 
galaxy surveys, such as SP SS (,Einasto et al...2003.) and 
APM ijEinasto et al.ll200^ . as wel l as individual struc- 
tures including the Ursa Major ("KoDvlova & KodvIov) 
Shapley l)Quintana,_^Carrasc o . fc Rciscncgger 
20001) Corona Bor'ealiris"^^ etTall Il998l) . Aquar- 
ius iCaretta et al l I2002D . and North Echptic Pole 
ijMullis et alJl2001(l superclusters. 

4. DISCUSSION 

We have obtained the largest spectroscopic sample to 
date for a supercluster at z ^ 0.9. The initial results 
suggest that the four component clusters are part of a 
single, filamentary structure, with the individual clus- 
ters being massive (Abell Richness class 0-2). As seen 
in Figure O we may not yet have reached the outer lim- 
its of this structure, and future wide-area imaging and 
additional spectroscopy will be required to measure its 
full extent. Already our dataset is comparable in size to 
studies of local superclusters (see §3.2). This structure is 



by far the best studied of the very few known superclus- 
ters at z > 0.3 llConnollv et alJll996t iRosati et al1ll999t 
iCavazzi et alJl2004D . 

The spectroscopic data, especially from DEIMOS, are 
of sufficient resolution and quality to measure equiva- 
lent widths for both emission and absorption features. A 
large fraction of the galaxies in the cluster show evidence 
for star- formation (most notably strong [Oil] emission). 
Examining the correlation of such features with galaxy 
location will yield insights into the dynamics and his- 
tory of the supercluster system and shed light on the 
role of cluster-scale physical processes (e.g., shocks and 
ram-pressure stripping) in galaxy star formation. Pro- 
posed HST imaging and Chandra X-ray mapping of this 
structure will allow us to compare galaxy morphological 
and spectroscopic properties and correlate optical sub- 
structure with the location and temperature of the hot 
gas component. These data, combined with numerical 
simulations, may be able to predict future mergers and 
perhaps even the end-state of this system. This survey 
will provide a view of a young supercluster rivaling the 
detailed studies of local structures. 

We dedicate this letter to the memory of Bev Oke, 
an extraordinary scientist, collaborator, and friend. We 
thank the anonymous referee, C. Fassnacht and G. 
Squires for helpful comments and contributions. Data 
presented herein were obtained at the W.M. Keck Ob- 
servatory, which is operated as a scientific partnership 
among the California Institute of Technology, the Univer- 
sity of California and the National Aeronautics and Space 
Administration. The Observatory was made possible by 
the generous financial support of the W.M. Keck Foun- 
dation. This research was supported in part by STScI 
grant HST-GO-08560.05-A. 
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